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Abstract 

We study the color-dependent confining forces between two quarks by the quenched lattice 
simulations of Coulomb gauge QCD. The color-singlet and color-antitriplet instantaneous potentials 
yield attractive forces. The ratio of the string tensions obtained from them is approximately 2 and 
have little volume dependence. Meanwhile, the color-octet and color-sextet channels give a minor 
contribution for two-quark system. We finally find that the infrared self-energy of the color- 
nonsinglet channels diverges in the infinite volume limit; however, the degree of the divergence on 
the finite lattice can be understood in terms of color factors. 
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I. INTRODUCTION 



The long-distance color-dependent forces among quarks and gluons are the key quanti- 
ties for the understanding of the internal structure of hadron as well as color confinement 
dynamics. The quarks and gluons are described in the quantum chromodynamics (QCD), 
and have a color charge based on the SU(3) group, with which the quark combination yields 
many color-dependent forces including both attractive and repulsive forces. They make the 
hadron structure (quark bound state) more complicated. Therefore, it is essential to know 
the basic behavior of the color-dependent force at short and large distances. 

The color-dependent force will be important when one investigates the multiquark hadron 
(made of more than four quarks) and the exotic meson, etc. Although many candidates of 
those particles have been recently reported by experimental groups y] , it is still a hard 
task to theoretically understand those new particles. They may consist of many quarks 
with various color-dependent forces that make their internal structure more complicated 
than the existing baryons. The quark interaction at short distances can be characterized 
by a coefficient (A ■ A) of a one-gluon exchange potential. This is a basic assumption if we 
construct a quark model to deal with the baryon system. However, owing to the quark 
confinement, it is not obvious how quarks behave for large quark separation. We thus need 
a lattice simulation to obtain a nonperturbative feature for the long-range color-dependent 
force. 

However, the lattice study along this line is sparse, although there are many lattice 
studies about the qq potential obtained by the gauge-invariant Wilson-loop operator. This 
operator mixes color-singlet with color-octet contributions; therefore, we cannot separately 
extract the color-octet potential from the Wilson loop. In spite of this, the color-octet 
process has been discussed in the detailed analysis of the J/ip photoproduction to explain 
the deviation between the experimental data and the theoretical predictions . On the 

other hand, there is no lattice calculation for qq sector as well; however, the diquark (color-3* 
qq) is also an important ingredient on hadron physics. The existence of the diquark has been 
believed by many physicists and applied to work out unsolved problems 0, ^, I^, 1^, 11, 12|. 
Note that the qq potential in the lattice theory is not formulated gauge-invariantly. 

In this study, we employ a color-dependent Polyakov line correlator (PLC) with the 
Coulomb gauge fixing in order to clarify the color-dependent forces 13|, llJ| . In the quark- 
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gluon plasma phase, there have been already some numerical studies of the color-screened 



color-dependent force with the Coulomb and Landau gauge fixings 
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16|, 



tigating the confining color-dependent forces we use Coulomb gauge QCD 



background of which has been recently studied well 
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17| 
l3 



18| . Inves- 



theoretical 



2l|. In this theory, the PLC 



potential can be separated into the vacuum polarization part and the color-Coulomb in- 



stantaneous part; the 



confinement phase 
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after relatively produces us a clear numerical signal even in the 
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251]. The color-Coulomb instantaneous potential is not the 



same as the Wilson-loop potential. However, it yields an upper bound of a linearly-rising 
quark potential {26 1 and is the most important quantity on the Coulomb gauge confinement 
scenario |l9]. The infrared singularity of this potential is caused from accumulation of the 
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Faddeev-Popov (FP) ghost eigenvalue at the vanishing momentum 

Coulomb gauge QCD may provide us the most favorable framework when one considers 
the hadron phenomenological studies such as a constituent quark model. The Coulomb 
gauge as a physical gauge gives a positive-definite Fock space that is suitable for constructing 
effective Hamiltonians. Moreover, the color-Coulomb instantaneous potential without the 
vacuum polarization (or a retarded effect) is required to make a quark bound state in analogy 
with the quantum electrodynamics (QED). This formulation is applicable in the heavy-quark 
system with no retarded effect as well as in the light-quark system with the constituent 



quarks due to chiral-symmetry breaking 



29|. 



In this paper, we study the long-distance divergence behavior of the color-dependent 
forces. Because the color-nonsinglet state cannot exist in nature, singularities in such states 
will emerge on the lattice as a finite- volume effect. This has not been verified in the previous 
calculation [13], and thus we should investigate their divergence behavior on a variety of 
lattice sizes. We also discuss the color-factor dependence for the degree of the divergence 
and will confirm this point by the numerical lattice simulation. In section II, we briefly 
summarize Coulomb gauge QCD and the definition of the color-dependent potentials with 
the PLCs. In section III, we give lattice numerical results and fitting analyses. Section IV 
is devoted to our summary. 
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II. COLOR-DEPENDENT POTENTIALS 



A. Instantaneous potential in Coulomb gauge QCD 

Coulomb gauge QCD has been quantized through the Faddeev-Popov technique |3] 
and renormahzabihty of this theory has been also proved in terms of the Hamiltonian and 



Lagrangian formalism j20|, l2l|| . The use of the Coulomb gauge as a physical gauge leads us 
to classify transverse gluon modes and an instantaneous interaction, which is required to 
make quark bound states in analogy with QED. 

The Hamiltonian of QCD in the Coulomb gauge can be given by 



2 J V » V ' ■ * V // ■ 2 _ 
where E^"^, Bi and p are the transverse electric field, the transverse magnetic field and 
the color-charge density, respectively. The function V in the second term is made by the 
Faddeev-Popov (FP) operator in the spatial direction, M = —Dd = —(9^ + gA x d), 



V{x,y)= / £z 



(2) 



_M{x,z)^ '^^^M{z,y) 
From the partition function with the Hamiltonian Eq. ([1]), one can evaluate the time-time 
gluon propagator composed of the following two parts: 

g'{Ao{x)Ao{y)) = g^D,,{x - y) = V{x - y) + P{x - y), (3) 

where 

V{x-y)=g^{V{x,y))6{x,-y,). (4) 

Equation (jl]) is the instantaneous color-Coulomb potential at equal time and causes antis- 
creening, so that this potential should be a confining potential to attract quarks in hadrons 
and is the most important quantity on the Coulomb gauge confinement scenario. Note that 
Eq. (jl]) in the case of quantum electrodynamics (QED) as a nonconfining theory is identified 
as a Coulomb propagator {—1/df) or a Coulomb potential 1/r. 
Simultaneously, Eq. ([3]) stands for a vacuum polarization, 

P{x-y) = -9'{j V{x,z)p{z,x,)£z j V{y,z')p{z\y,)d'z'), (5) 

which brings about color-screening effect owing to the negative sign, producing the reduction 
of a color- confining force. Note that this term is associated with a quark-pair creation from 



vacuum when dynamical quarks exist. Consequently, we can consider that the instantaneous 
interaction is somewhat more classical than the time-dependent vacuum part. 



B. Color-Coulomb instantaneous potential 

One can define the color- dependent potentials on a lattice with the Polyakov line (PL) 
correlators In this study we separate the original potential obtained from the PL cor- 
relators into the color-Coulomb instantaneous and color vacuum polarization (retarded) 
parts with the Coulomb gauge fixing. Moreover, the color-Coulomb potential defined by 
the link-link correlator as will be described below gives clear signals into practical numerical 
calculation even in the quenched lattice simulations. 



We employ a partial-length Polyakov line (PPL) which can be defined as 22l. |23| 



L(f,T) = J]f/o(f,t), T = l,2,--- ,iVi. (6) 

t=i 

Here Uo{x,t) = exp{iagAQ{x,t)) is an SU{3) link variable in the temporal direction and a, 
g, Ao(x, t) and Nf represent the lattice cutoff, the gauge coupling, the time component of a 
gauge potential and the temporal-lattice size. A PPL correlator in the color-singlet channel 
is given by 

G^{R,T) = ^{Tt[L{R,T)L\0,T)]) , (7) 

where R stands for From Eq. ([7]) one evaluates a color-singlet potential on a lattice, 

Gi(R,T) 



For the smallest temporal-lattice extension, i.e., T = 0, we define 



(8) 



V{R,0) = -\og[G,iR,l)]. (9) 

Here V{R, 0) in the Coulomb gauge is assumed to be the color-Coulomb instantaneous 
potential Vcoui{R)- The V{R,T) in the limit T (yo becomes the usual Polyakov line 
correlator. These two potentials are expected to satisfy Zwanziger's inequality, Vphys{R) < 
Vcoui{R) {3], where Vphys{R) is the physical potential extracted from the Wilson loop. 
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C. Color-dependent potentials on a lattice 



We apply t 



two quarks 



le above discussion to the other SU{3) color-dependent potentials between 



14( 1 ■ A color-octet correlator on qq is given by 



Gs{R,T) = ^(TrL(i?,T)TrLt(0,T)>-^(TrL(i?,T)Lt(0,T)>, (10) 

and qq correlators in the symmetric-sextet and antisymmetric-triplet channels (3(8)3 = 6 ©3) 
are also given as 

Ge{R, T) = 1 {TtL{R, T)TrL(0, T)) + ^ {TtL{R, T)L{0, T)) , (11) 

Gs{R, T) = ^ {TrL{R, T)TrL(0, T)) - ^ {TrL{R, T)L(0, T)) . (12) 

In the same way as described in Eqs. ([8]) and ([9]) we obtain the color-dependent potentials 
in each color channel. 

The above four potentials are classified in terms of the color (the quadratic Casimir) 
factor on color SU{3) group in the fundamental representation: 

cL- = --, C- = -> ci = --, C',, = -, (13) 

qq 2 5 qq q1 qq g > qq ^i V / 

for color- singlet, color-octet, color-triplet and color-sextet channels, respectively. These 
coefficients appear as the proportional constant of the one-gluon exchange potential ^ 



D. Gauge fixing 

Since the color-dependent potentials defined by the PPL correlators are not gauge invari- 
ant, we have to fix a gauge. We use the Coulomb gauge realized on a lattice as 

3 



Max J2 I^eTr?// (f , t) , (14) 



by repeating the following gauge rotations: 



Ui{x,t) U^{x,t) = uj\x,t)Ui{x,t)uj{x + t,t), (15) 



^ Here we intend to discuss an irreducible representation for qq and qq sectors in relation to a so-called (A • A) 
quark model. Therefore, in this study, we do not calculate the gauge- invariant Wilson-loop potentials in 



higher representations for qq sector as introduced in Refs. [30l l3lU32ll33| 
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where uo G S'?7(3) is a gauge rotation matrix ^ and Ui{x,t) are spatial lattice link variables. 



34|. 



The thermalized lattice configuration can be gauge fixed iteratively 

Because the Coulomb gauge fixing does not fully fix a gauge one can still perform a 
time-dependent gauge rotation on the Coulomb-gauge fixed links, 

(16) 

Uo{x,t) uj^{t)Uoix,t)uj{t + l). 

Thereby, TrLTrL^ and TiLL constructed by PPL are not invariant under this transformation 
Eq. (fT6|) . Accordingly, when performing numerical simulations for the octet and two qq 
correlators with TrLTrL^ and TtLL, we should additionally implement a global temporal- 
gauge fixing on the Coulomb-gauge fixed links as 

Max ^ ReTrt/ct(f , t) under Eq. 1^, (17) 

x,t 

where V = N^NyNz is a spatial lattice volume. Note that this gauge fixing does not affect 
an intrinsic Coulomb gauge feature. 



E. Infrared divergence of the color-dependent potential 

In Coulomb gauge QCD, we find that infrared divergences cancel only for color-singlet 
interactions. Let us write the color-charge density for two-quark system as pa ~ T^5{x — 
xq) +T2S{x — yo), where are the generators of the color-5'?7(3) group. The instantaneous 
potential with R = xq — yo (Eq. ([2])) is given by 

V(R) = TfTl / Ji^^fiP}fMe'^-S, (18) 

where d{p) is the expectation value of the Faddeev- Popov operator and f{p) means the 
deviations of the factorization, e.g., f{p) = 1 if there is no quantum correction; this argument 



has been done in Refs. 



35 



36[ |. Here TfTg^ is reduced to each color factor that appeared in 
Eq. (fT3|) . On the other hand, the self-energy interaction not depending on the distance R 
is written as 



^ Here we used uj — e'"^»"^', where the parameter a is chosen suitably depending on the lattice size, etc. 
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where the Casimir invariant (T^)'^ = 4/3 in the fundamental representation of the SU{3) 
group. 

The infrared (not ultraviolet) divergence emerges in both Eps. 0181) and (I19p . If the 
term d'^{p)f{p) behaves as (1/ ^/py • 1/p that would be responsible for the linear confinement 



V ~ as has been analyzed in the Refs. 36|, 1371], then the infrared singularities unrelated 
to the linear potential arise from two terms: 

V^'iR) = 4vr(TfT2^) / dp-, S^^ = AniTn' / dp-. (20) 
Jo P Jo P 

Consequently, they are completely cancelled in the case of the color-singlet representation 

since (TfTg) + (T")^ = (—4/3) +4/3 = 0. Meanwhile, the other cases are proportional to the 

following factors: 3/2, 2/3, 5/3 for 8, 3* and 6, respectively, implying that the color-sextet 

channel may diverge most strongly. 



III. NUMERICAL RESULTS 



A. Simulation parameters and statistics 

We carry out SU{3) lattice gauge simulations in the quenched approximation to calculate 
the color-decomposed PPL correlators. The lattice update is done by the heat-bath Monte 
Carlo algorithm with a plaquette Wilson gauge action. The lattice configuration numbers 
for the 18^, 24^ and 32^ lattices are 600, 700 and 320; additionally, in order to investigate the 
volume dependence, we added the 8^ and 12'* lattices with 200 configurations. The lattice 
coupling constant (3 for all the lattices is fixed to 5.9 corresponding to the lattice cutoff 
a ~ 0.12/m jssl. 



B. Color-dependent potentials 

Figure [1] shows numerical results for the color-Coulomb instantaneous potential V{R,0), 
in the color- singlet, color-octet, color-sextet, color-triplet (antisymmetric) channels. We find 
that both the color-singlet Vi and color-antitriplet V3* yield attractions at all distances, and 
in particular, are linearly rising potential at large distances. On the other hand, the color- 
octet Vs and color-sextet Vq potentials are repulsive forces although the variation of those 
potentials on distances is small. 
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p=5.9,a=0.123 




FIG. 1: Lattice numerical results of color-Coulomb instantaneous potentials between two quarks. 
{P = 5.9, a ~ 0.12/m) 




FIG. 2: Effective string tensions for the color-singlet and color-antitriplet channels in lattice units. 
Open (Filled) symbols are the color-singlet (antitriplet) string tensions. Here the error bars are 
estimated in the error propagation. The effective string tensions at large R are strongly affected by 
the periodic boundary condition of a finite lattice; in contrast, owing to the increase of the lattice 
size, they become more stable even at large R. 
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In order to investigate the magnitude of the string tensions, we calculate an effective 
string tension for the color-singlet and color-antitriplet channels. The confining potential 
can be described by the function V{R) = C + KR + A/R where C is a self-energy constant 
term, K corresponds to the string tension, and the last term is the Coulomb term. Here the 
effective string tension is defined as K = V{R + 1) — V{R) in lattice units, which should be a 
constant for large quark separations if it is a confining potential with a finite string tension. 
In Fig. [2] we find that the Ks become stable over approximately i? = 3 ~ 4 as the lattice size 
increases. When we use the data for i? = 3 — 6 (i? = 3 — 5 for the 18^ lattice), we plotted in 
Fig. [3] the ratio K1/K3*, which is found to be close to C1/C3* = 2; for example, we obtain 
the Ki = 0.218(2) and K^, = 0.127(13) from the 24^ lattice. Note that our definition of the 
color-Coulomb instantaneous part on a lattice in terms of the PPL correlator (Eq. (^) does 
not completely exclude a vacuum polarization effect. ^ 

For the color-octet and sextet channels it is hard to discuss a precise value of both string 
tensions under the present statistical accuracy. In practice, from 24^ lattice data and under 
the fitting range = 3 — 6, we obtain the octet string tension Kg = 0.0014(222) and the 
sextet string tension Kq = —0.018(19) although the volume dependence does not appear 
large as shown in Fig. 4. This indicates that the long-range behavior of the octet and sextet 
channels is quantitatively a minor contribution to two-quark system. However, in order to 
work out more quantitative measurement for the repulsive channels (and also the attractive 
channels), one may need enough statistics on a larger lattice. 



C. Divergence of color- no nsinglet potentials 

Here we consider the infrared divergence property of the color-dependent potentials while 
the string tensions obtained from them have little volume dependence as shown in the pre- 
vious section. The short- distance Coulomb term ~ 1/i?, which is not related to the color 
confinement, may not matter on this argument. Figure [5] shows our nonperturbative numer- 
ical results, by which we find that the color-singlet potential has little volume dependence 
and the absolute value of the color-nonsinglet potentials increases with the lattice volume; 
this tendency remains at the distances R = 3 and 6 as displayed in Fig. [6l Thereby we 
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A possible way to improve this discussion may be to construct directly the instantaneous potential by the 



Faddeev-Popov propagator that has an infrared singularity 271 128| 
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FIG. 3: Vertical axis stands for ratios Ki/K^* of the string tensions in the singlet and antitriplet 
channels and horizontal axis is the spatial lattice size. 




FIG. 4: Effective string tensions for the color-octet and color-sextet channels in lattice units. 

conclude that the color-nonsinglet potentials diverge regardless of the distance. 

Moreover, in order to estimate the degree of divergence for the color-nonsinglet channels, 
we assume the fitting function f{R) — C -\- KR for the potential data in the distance 
it! = 3 — 6 (long range). The corresponding self-energy term Cs fitted well are plotted in 
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Fig. [71 this graph shows that the divergence contribution increases monotonically with the 
volume. Therefore we adopt the hnear function y = Dx + d where D may correspond to the 
constant meaning the degree of divergence. The fitted results are summarized in Table I, 
which gives the conclusion that the magnitude of the divergence on the finite lattice depends 
on the ratio being 4 : 9 : 10 for 3*, 8 and 6 as observed in the section II. E. However, note 
that they will diverge equally to infinity in the large volume limit, as shown in Figs. [5] and 
[71 This indicates that the color-nonsinglet quarks cannot exist independently; finally, they 
will become a color-singlet state with finite energy to compensate insufficient color degrees 
of freedom. 

This numerical result does not contradict the dual-Ginzburg-Landau picture, in which 
one understands that the color-singlet flux between two quarks is shrunk like a string while 
the color-nonsinglet flux will radiate rather than produce a closed string. A singularity of the 
color flux not making the closed string can be expected to emerge as the volume dependence 
on the lattice gauge theory. 




FIG. 5: Volume dependence of the color-dependent instantaneous potentials at the distance R = 1 
scaled by the V^{1, 0) on the 8^ lattice. The lattice sizes (L) used here are 8, 12, 18, 24 and 32. The 
variation of the singlet potential with the lattice volume is little while the other color-nonsinglet 
potentials diverge in the infinite volume limit. 
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FIG. 6: Volume dependence of the color-dependent instantaneous potentials (/3 = 5.9) at the 
distances R = 1, 3, 6. 
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FIG. 7: Volume dependence of the self-energy term for color-nonsinglet potentials. The self-energy 
term diverges in the infinite volume limit. 

IV. SUMMARY 

We have studied the long-distance color-dependent forces between two quarks in the 
quenched S'[/(3) lattice simulation with the Polyakov line correlator. We here focus the 
color-Coulomb instantaneous term in Coulomb gauge QCD, which has been discussed in 
the Gribov-Zwanziger confinement scenario and is required to make the hadron bound state 
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TABLE I: Table of fitted results to estimate the degree of divergence in color-nonsinglet channels. 
y = Dx + d is assumed as a fitting function to fit the data of Fig. [71 





3* 


8 


6 


D 


0.044(06) 


0.092(14) 


0.098(14) 


d 


0.42(07) 


0.38(10) 


0.39(10) 


X^/ndf 


0.19 


4.20 


4.63 



consisting of quarks. 

Our numerical simulation shows that the color-singlet qq channel as well as the color- 
antitriplet qq (diquark) channel causes a linearly confining potential at large distances. The 
other color-octet and color-sextet channels at large distances yield also repulsive forces. In 
addition, we find that the string tensions in the color-singlet and color- ant itriplet channels 
have not significant volume dependence 

We also investigated the infrared divergence of color-nonsinglet potentials and find that 
the divergence on the finite lattice seems to be proportional to the color (Casimir) factor; 
however, the divergence of the color-singlet channel is not left over. In the infinite volume 
limit the color-dependent potentials except the color-singlet channel will diverge; in particu- 
lar, the color-sextet channel diverges most strongly as expected. This conclusion is consistent 
with the dual-Ginzburg-landau picture of the color confinement. 

This approach, in which we take notice of the color-Coulomb instantaneous part, may 
be suitable for further work to investigate three-quark color-dependent forces for the un- 
derstanding of baryons as well as new multiquark particles. It is also necessary to study 
how the vacuum polarization term affects the instantaneous forces because it is reported in 
Refs. 



391] that the octet (adjoint) channel calculated by the full-length Polyakov line 



correlator with the vacuum polarization gives the complicated distance dependence at large 
distances. 
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VI. APPENDIX 

In Eq. f fTSj) . if the Coulomb kernel K{p) produces a confining force, it is reported in Ref. 
371] that the ghost factor d{p) and the Laplacian factor f{p) behave as approximately l/i/p 
and 1/p, respectively. Thus we assume that K{p) is 1/p^. Omitting T^T^ term one also 
calculates Eq. 



^ ^ ' 27r p2 

2tt pI 







/ dcosO r p^dpO^e'P^'"'' 
J-i Jo P 



^ r dp^ie'P'' - e-'P"") (21) 
iR Jq p 

27C r ,^K{p)^^^^ 



dp^^e' 

iR J-OO P 

= % I dp^{l + ipR + hipRf + \{ipRf + ■■■)• 
iRJ-oo P 2 3 

We here do the Taylor expansion in the infrared limit p = for e*^^. In the above equation, 
the first term disappears because it is a odd function on p; the second term proportional to 
l/p"^ causes the infrared divergence; the third term produces a (linear) potential; the other 
terms are irrelevant. Finally, the second term relating to the divergence is 

POO -r 

V^'iR) = 2n{T,^T^) / dp- 

■^788 ^ (22) 

= 47r(Tfr|) / dp-. 

Jo P 



In the same way, the relevant contribution from the self-energy of Eq. (1191) is also given by 

POO 1 

E^^ = 47i{Ttf / dp-. (23) 
Jo P 
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